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ABSTRACT 

Observations of transversal coronal loop oscillations very often show the ex- 
citation and damping of oscillations in groups of coronal loops rather than in 
individual and isolated structures. We present results on the oscillatory proper- 
ties (periods, damping rates, and spatial distribution of perturbations) for reso- 
nantly damped oscillations in a system of two inhomogeneous coronal slabs and 
compare them to the properties found in single slab loop models. A system of 
two identical coronal loops is modeled, in Cartesian geometry, as being composed 
by two density enhancements. The linear magnetohydrodynamic (MHD) wave 
equations for oblique propagation of waves are solved and the damping of the 
different solutions, due to the transversal inhomogeneity of the density profile, 
is computed. The physics of the obtained results is analyzed by an examination 
of the perturbed physical variables. We find that, due to the interaction be- 
tween the loops, the normal modes of oscillation present in a single slab split into 
symmetric and antisymmetric oscillations when a system of two identical slabs 
is considered. The frequencies of these solutions may differ from the single slab 
results when the distance between the loops is of the order of a few slab widths. 
Oblique propagation of waves weakens this interaction, since solutions become 
more confined to the edges of the slabs. The damping is strong for surface-like 
oscillations, while sausage body-like solutions are unaffected. For some solu- 
tions, and small slab separations, the damping in a system of two loops differs 
substantially from the damping of a single loop. 
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INTRODUCTION 



The solar corona, the outermost part of the solar atmosphere, is able to support a va- 
riety of MHD waves in its different magnetic and plasma configurations. In the last years, 
particular attention has been devoted to the phenomenon of trans yersal coronal loop oscil- 



lations, first observed by instruments on-board TRACE spacecraft (lAschwanden et al.l Il999 



Nakariakov et al. 1999 : Aschwanden et al. 2002 : Schrijver et al. 20021 ). These oscillations, 
triggered mainly by nearby flares or filament eruptions, consist of lateral displacements of 
the loops while their foot-points are fixed to the denser photosphere. The characteristic 
periods are of the order of minutes and oscillations are quickly damped in a few periods. 



From a theoretical point of view, these oscillations have been interpreted by lNakariakov et al. 



( 119991 ) as the fundamental fast MHD kink mode of a flux tube. This solution involves a 
global motion of the whole structure that displaces the axis of the tube. The identification 
of observed oscillations with th e oretic al MHD wave so l utions is the key for coronal seismol- 
ogy, first suggested by lUchidal (Il970l ) ; [Roberts et al.l (119841 ) , and has allowed preliminary 
estima tions of the magnetic f i eld strength, transport coe f ficients, and plasrn a-/? in coronal 
loops ( Nakariakov et al.lll999l : iNakariakov fc Ofmanll200ll : lAschwanden 12005). Some recent 
examp l es of the application o f coronal seisrnology c an be found in INakariakov fc Ofman 
(l200ll):lGoossens et all J2002h : lAndries et all J2005h : IVerwichte et al.l J2006h : lArregui et al. 



( l2007al ). As for the nature of the damping mechanism(s) there is little consensus yet and 



several mechanisms are currently u nder study, such as no n-ideal effects, lateral wave leakage 
due to the curvature of the loo ps (iTerradas et al.l l2006l) . mechanisms based on the topol- 
ogy of the magnetic field lines (ISchrijver fc Brownl |2000| ) . and resonant co nversion of wave 
energy, due to the non-uniformity of the magnetic-pla s ma co nfiguration (iHoUweg fc Yang 



19881 : iRuderman fc Robert] I2OO2I : ICoossens etaPbood bood). Foot-poin t leakage through 



the chromospheric density gradient (IDe Pontieu et al 



2001 



OfmanI I2OO2I ) has been invoked 



for the damping of Alfven waves and could be ascribed to the damping of transverse oscil- 
lations, as well. 

Observations of transversal coronal loop oscillations very often show the excitation and 
damping of motions in groups of coronal loops rather than in single, isolate d structures (see 



for ex ample the event recorded by TRACE on 2001 April 15, analyzed by IVerwichte et al. 



2004 ). However, most of the theoretical models used for their study are based on single 
loop models. The interaction between the motions of adjacent coronal loops may change the 
oscillatory properties (periods, damping times, and spatial distribution of perturbations) of 
the observed event. For this reason, it is important to fully understand the theoretical prop- 
erties of MHD waves in multiple-loop structures. The interest of these studies is two-fold. 
On one hand, to see under which conditions these differences in the oscillatory properties 
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are important and, on the other hand, if important, to make rehab le predictions and deter- 
minations of unknow n physical parameters in the corona. For example, when analyzing the 
2001 April 15 event, IVerwichte et al.l (120041 ) report on damping times that are longer than 
expected from previous observations of transverse oscillations in isolated loops. The fact 
that the studied coronal loops belong to a group could be behind this result. 

Several authors have previ o usly investigated the properties of MHD waves in multi- 
structures. iBerton fc Hevvaertsi (119871) stud i ed th e pro perties of the normal modes of oscil- 



lation of a periodic medium. iBogdan fc Foxl (Il99ll ) and lKeppens et al.l (Il994j ) considered the 
scattering and absorp t ion o f aco ustic waves by bundle s of rn agnetic flux tubes with sunspot 
properties. iMurawskil (Il993l ) and lMurawski fc Roberts! (119941 ) studied, by means of numerical 
simulations, the propagation of fast waves in a system of two slabs unbounded in the longitu- 
dinal direction. The coll ective nature and properties of oscillations in multi-fibril Cartesian 
systems was ana. l yzed b y lDiaz et al.l (120051 ). in the context of prominence oscillations. More 
recently, lOfman 



J2005h perfor med numerical simulations f or the damping of a bundl e of 



four loops, iMarcu et al.l (120061 ) have extended the study by iBerton fc Hevvaertsi (119871) by 
including a steady motion in the equilibrium and lOgrodowczyk fc Murawskil (l2007l ) have 
performed numerical simulations of impulsively exc ited magneto-sonic waves in two parallel 



solar coronal slabs. Also recently, iLuna et al.l (120061 ) have explored the normal modes and the 



time-dependent evolution of fast kink MHD perturbations in a system of two coronal loops, 
modeled as plasma slabs. They find a splitting of modes into symmetric and antisymmetric 
kink-like oscillations. The interaction between two coronal loops produces a variation on 
the frequency and spatial structure of the eigenfunctions, that depend now on the distance 
between the loops. 

In this paper, we present the oscillatory properties of fast MHD waves in a system of 
two Cartesian slabs when both oblique propagation of perturbations and transversal inhomo- 
geneity of the medium are included. These two ingredients produce the resonant damping of 
fast waves. A particu lar study of the oscilla tory properties of a similar physical system was 
recently presented by lArregui et al.l (l2007bl ). who considered two very close slabs to analyze 
the effect of the possibly unresolved internal structuring of coronal loops on their transverse 
damped oscillations. Here a more in-depth analysis is performed. The reso nant damping 



of fast MHD oscillations in a single Cartesian slab has been described by lArregui et al. 



(l2007cl ). In this paper, a second density enhancement is in cluded and a n ana logous analysis 
is presented. This work, thus, constitutes an extension to iLuna et al.l (120061 ) to the case of 



oblique propagation and damping by resonant absorption and an extension to lArregui et al. 



(l2007cl ) to a system of two interacting slabs. Although the model adopted in this work is 
an oversimplification, and is far from being an accurate representation of real coronal loop 
systems, this simplicity allows a detailed study of the problem. 
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The layout of the paper is as follows. In § |2l the equilibrium configuration, the linear 
MHD wave equations for coupled fast and Alfven waves, and a description of the numerical 
method are presented. Then, §[3]describes the results of our study. First, in § 13.11 the normal 
mode properties of undamped fast waves for oblique propagation are described. Then, in 
§ 13. 2[ the damping properties of the solutions and the physical interpretation of the results 
are discussed. Finally, in § HI our conclusions are drawn. 



2. EQUILIBRIUM MODEL, LINEAR MHD WAVE EQUATIONS, AND 

NUMERICAL METHOD 



We model the equilibrium magnetic and plasma configuration of a system of two coronal 
loops by means of a one-dimensional model in Cartesian geometry. The magnetic field is 
straight and pointing in the z-direction, B = i? e^. For applications to the solar corona, 
it is a good approximation to consider that the magnetic pressure dominates over the gas 
pressure. This classic zero plasma-/? limit implies that the magnetic field is uniform and that 
the density, p{x), or Alfven speed, va{x), profiles can be chosen arbitrarily. For simplicity, 
we consider a system of two identical coronal loops. This system is then modeled by defining 
a particular equilibrium density profile in the x-direction (see Figured]), with two density 
enhancements of half-width a located at ixg. The density in each of the slabs is uniform. 
Pi, and connected to the uniform coronal environment, with density Pe, by transitional non- 
uniform layers of thickness /. The explicit expression for the considered equilibrium density 
is (for a; > 0) 
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with p{—x) = p{x). The particular dependence of the density profiles at the non-uniform 
transitional lay ers has no qualitative i mport ance on the obtained results . We have chosen 
them, following iRuderman fc Roberts! (120021 ) and IVan Doorsselaere et al.l (120041 ). as 
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and 



This equilibrium configuration is the same as the one used by lArregui et aLl (j2007bl ) 
for the study of the influence of the internal density structuring of a coronal loop on the 
properties of their damped oscillations. Their analysis is based on an equilibrium model in 
which the two slabs are very close to each other and contains a comparison of the results 
to the ones obtained with similar s ingle slab models. This equilibrium model is also an 
extension of the one considered by iLuna et al.l (120061 ) . with the addition of non-uniform 
density layers at the edges of each slab. 

In order to study small amplitude oscillations of the previous equilibrium, the linear 
resistive MHD equations, with constant magnetic diffusivity, 77, are considered. As the equi- 
librium configuration only depends on the x-direction, a spatial and temporal dependence 
of the form Qxp^^^^~''y^~^''^^ is assumed for all perturbed quantities, with cu = ujji + loji the 
complex frequency and ky and the perpendicular and parallel wavenumbers. For reso- 
nantly damped solutions, the real part of the frequency gives the period of the oscillation, 
P = 2-71 /ujfi, while the imaginary part is related to the damping time, Td = ^/uji. The pho- 
tospheric line-tying effect is then included by selecting the appropriate parallel wavenumber. 
This leads to the following set of ordinary differential equations for the two components of 
the velocity perturbation, and Vy, and the three components of the perturbed magnetic 
field, bx, by, and bz, 
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These equations, together with the appropriate boundary conditions, define an eigenvalue 
problem for fast and Alfven normal modes. As the plasma-/3=0, the slow mode is absent 
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and there are no motions parallel to the equilibrium magnetic field, Vz = 0. In this paper, 
solutions to these equations are obtained by performing a normal mode analysis. When 
only perpendicular propagation of perturbations is considered {ky 7^ 0, / = 0), the solutions 
represent undamped fast MHD waves with oblique propagation, in addition to Alfven waves. 
If non-uniform transitional layers are included (/ 7^ 0), the solutions are damped. In this 
case, solutions are in general difficult to ob tain analytically. Numerical approximations are 
then obtained using PDE2D (ISewelll l2005l ) . a general-purpose partial differential equation 
solver. The code uses finite elements and allows the use of a non-uniformly distributed grid, 
which is needed in order to properly resolve the large gradients that arise in the vicinity of 
resonant layers. This code has been used successfully i n previous studie s that i nvolved similar 



comput ations of the damping rate by, for example, iTerradas et al.l (120061 ): lArregui et al. 



f l2007bl l3l As for the boundary conditions, we impose the vanishing of the perturbed velocity 
far away from the two-slab system, hence v as a; — i> ±cxd. 



ANALYSIS AND RESULTS 



Solutions to 
been ob tained by 



quations p])-(171) for the simplest case, when ky = Q and / = 0, have 



Luna et al.l (120061 ). Fast normal modes are then decoupled from Alfven 

modes. Luna et al. ( 20061 ) find that the system is able to support a variety of oscillatory 



modes. In particular, they concentrate on the splitting of the kink transversal oscillation 
of a single slab, giving rise to two solutions, symmetric and antisymmetric with respect to 
X = 0, in which the two slabs oscillate in the kink mode in phase or anti-phase with the same 
frequency. The fundamental mode is symmetric and its frequency is constrained between 
two limiting values. In the limit of no separation between slabs, it is equal to the kink 
mode frequency of a single slab with double width. As the distance between the slabs is 
increased, the interaction weakens and, for large separations, the two slabs do not "feel" 
each other and the frequency of the whole system is the same as the kink mode frequency 
of a single slab. The range of separation between slabs for which the interaction is strong is 
mainly determined by the parallel wavenumber, fc^, which determines the confinement and 
drop-off rate of the eigenfunctions in the coronal medium. These authors also find that the 
antisymmetric mode is not present for all the values of the distance between the slabs and 
becomes leaky when this distance is smaller than a given value, that depends on the adopted 
value for the parallel wavenumber. 



When compared to iLuna et al.l (120061 ) . there are two new ingredients in the physical 
model under consideration in this paper. On one hand, the non-uniformity of the equilib- 
rium density at the edges of the slabs. On the other hand, the inclusion of perpendicular 
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propagation of perturbations. These two ingredients produce the resonant couphng of fast 
modes to Alfven modes and, as a result, the global modes of the system are damped. The 
damping of oscillations in a sys tem of coronal loops has been reported in observations ana- 
lyzed by IVerwichte et al.l (120041 ) and is likely to differ from the theoretical values obtained 
with isolated coronal loop models. 



In § 13. H we first describe the normal mode properties of oscillations of our two-slab 
system when oblique propagation of perturbations is included. Then, in § 13.21 the damping 
of oscillations by resonant absorption and the properties of the eigenfunctions are detailed. 



3.1. Normal Modes in a Two-Slab System with Oblique Propagation 



We first consider solutions to Equations Q-© including oblique propagation of waves 
{ky 0) and in the absence of non- uniform transitional layers (/ = 0). Also, rj = is set 
in these equations. As the equilibrium is piece-wise constant there is no resonant coupling 
to Alfven waves and solutions represent undamped oscillations. These solutions can be 
obtained by solving an analytical dispersion relation. lArregui et al.l (l2007d ) have shown, for 
a single Cartesian slab, that oblique propagation of waves affects the nature and properties 
of kink- and sausage-like solutions. The fundamental kink mode, which is body for values 
of ky below certain value, becomes surface beyond that particular value. There is a sausage 
surface wave, with its frequency always below the internal cut-off frequency, which only 
exists when ky 0. The phase speed of these two solutions goes asymptotically to the kink 
speed for quasi-perpendicular propagation {ky ^ k^). On the other hand, the sausage body 
solution, which is leaky in the long wave limit {k^a <^ 1), becomes trapped when a non-zero 
ky is included. This solution is always above the internal cut-off frequency and, hence, keeps 
its body character in the limit of quasi-perpendicular propagation, due to the existence of 
the sausage surface solu tion below the interna l cut-off frequency. All these features are well 
displayed in Figure 2 of lArregui et al.l (l2007d ). 



When a second density enhancement is included, the kink, sausage surface, and sausage 
body solutions, described by lArregui et al.l (l2007d ). split into symmetric and antisymmetric 
solutions with respect to x = 0, giving six solutions that we name symmetric and antisym- 
metric; kink (with the two slabs oscillating in the kink mode), sausage surface (with the two 
slabs oscillating in the sausage surface mode), and sausage body (with both slabs oscillating 
in the sausage body mode). 

When finding our solutions, we consider fixed values for the density contrast of the two- 
slab system, pj/pe=10, and for the longitudinal wavenumber, kza=n/50. For the observed 
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kink oscillations with a wavelength double the length of the loop, this corresponds to a 
ratio of length to width L/2a = 25. These are typical values for observed coronal loops 
for which transverse oscillations have been reported. We have studied the dependency of 
the frequency of these solutions as a function of the position of the slabs, related to the 
distance between the loops, and as a function of the perpendicular wavenumber. Figure [2] 
displays the obtained results. The splitting of the solutions can clearly be seen for the six 
modes. For the kink symmetric and antisymmetric solutions, the range of distances for which 
the interaction is strong and the deviation from the kink mode frequency of a single slab 
significant, is large in the absence of perpendicular propagation {ky = 0). Even for distances 
of the order of 100a, the two solutions do not converge to that frequency. When perpendicular 
propagation is included the single slab kink mode frequency is approached for much smaller 
distances between slabs, in such a way that only for distances of the order of a few times 
a the symmetric and antisymmetric solutions have frequencies that differ significantly from 
the kink frequency of a single slab. Also, the increase of the perpendicular wavenumber 
allows the existence of the kink antisymmetric mode for smaller distances between slabs, 
in comparison to the ky = case, when it becomes leaky below xq ~ 26a. As for the 
sausage surface and sausage body, symmetric and antisymmetric solutions, a similar result 
is obtained. They are plotted only for values of ky because the surface solutions only 
exist when ky and the body solutions are leaky below a certain value of ky. We can also 
appreciate the splitting of the single slab solutions and the dependency of the interaction 
between the slabs and, hence, the deviation from the single slab frequencies as a function 
of the slab separation and perpendicular wavenumber. It is interesting to note that the 
fundamental solution in the case of kink and sausage surface solutions is the symmetric one, 
while it is the antisymmetric one in the case of the sausage body solutions. The right panels 
of Figure [2] display the frequencies of the six solutions as a function of ky for a fixed value of 
the distance between the slabs. As we increase the value of the perpendicular wavenumber 
both symmetric and antisymmetric solutions converge to the same frequencies and approach 
those of the kink, sausage surface, and sausage body solutions of a single slab in the limit of 
quasi-perpendicular propagation. 

We have computed the difference between the frequency of the symmetric and anti- 
symmetric solutions of the two-slab system and the kink and sausage solutions of a single 
slab. For the kink solutions, and considering for example xo = 2a, the frequency of the 
symmetric mode is 15% smaller than the kink mode frequency of a single slab, when ky = 0. 
This difference decreases to 3% for kyU = 1. When the two-slab system is oscillating in the 
antisymmetric mode, its frequency is 13% larger that the kink mode frequency for kya = 0.5, 
while it is only 4% for kya = 1. As for the sausage surface solutions of the two-slab sys- 
tem, a deviation from the sausage surface frequency of a single slab around 5% is found, in 
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the case of the symmetric solution with kya = 0.3 and xq = 2a. Finally, for the sausage 
body solutions of the two-slab system with xq = 2a, the maximum difference occurs for the 
antisymmetric solution with kya = 0.5 and is around 10%. 

The results displayed in Figure [2] can easily be interpreted by analyzing the spatial 
distribution of the corresponding eigenf unctions. Figures |3l HI and [5] show the transversal 
component of the perturbed velocity and the total pressure perturbation for the kink, sausage 
surface, and sausage body, symmetric and antisymmetric solutions, for different values of 
the perpendicular wavenumber. In Figure [3] the solutions for the kink solutions are plotted 
for three values of the perpendicular wave number. For ky = 0, the transversal velocity is 
maximum at the internal part of the slabs, while it decreases exponentially outside. The 
total pressure perturbation has an antisymmetric profile in each slab, with a zero value at 
the center. When oblique propagation is included, there is an improved confinement of the 
eigenfunctions wi th a sharper drop-off rate in the external medium, a behavior also found by 



Diaz et al.l (120031 ). in the context of prominence fibril oscillations. Also, the character of the 
solutions inside the slabs change, and the two solutions become surface-hke solutions, with 
a decreasing amplitude of Vx inside the slabs as ky is increased. This improved confinement 
of the solutions weakens the interaction between the slabs and results in a frequency of the 
two-slab kink solution closer to that of the kink mode of a single slab. Regarding the sausage 
surface and body-like solutions (Figures [Hand [5]), we can appreciate that they have a similar 
spatial distribution of eigenfunctions. The main difference has to be found in the total 
pressure perturbation profiles, that are maximum at the center of each slab, for the sausage 
body solutions, while they peak at the edges of the slabs for the surface-like counterparts. 
An increase in the perpendicular wavenumber also produces an improved confinement, in 
the external regions, for these four types of modes. This also helps explaining why the 
interaction between the slabs weakens for increasing perpendicular wavenumber. However, in 
contrast to the kink solutions, now the transversal velocity component inside the slabs is not 
affected by the increase in ky and the solutions keep their surface or body character. Finally, 
oblique propagation of waves affects the magnitude of the total pressure perturbation of the 
sausage surface, symmetric and antisymmetric solutions (Figure H]), but not very much that 
of the sausage body, symmetric and antisymmetric, modes inside the density enhancements 
(Figure E]). 

In summary, the inclusion of perpendicular propagation of perturbations changes sig- 
nificantly the oscillatory properties of fast MHD waves in our two-slab system. First, it 
allows the existence of two sausage surface solutions that are not present when ky = 0. 
Also, it allows the existence of trapped sausage body solutions for slab separations for which 
they are leaky when ky = 0. Second, it produces an improved confinement of the solutions 
that, for the six modes, are better confined to the slabs, with a sharper drop-off rate in 
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the coronal medium. This produces changes in the frequencies of the six modes and varies 
significantly the range of slab separations for which the difference between the symmetric 
and antisymmetric frequencies and the single slab frequencies are important. In the case of 
the kink solutions, the increase of the perpendicular wavenumber also changes the character 
of the solutions inside the slabs, and the modes become surface-like. Therefore, for a fixed 
value of the parallel wavenumber, both the separation between the loops and the value of 
the perpendicular wavenumber determine the strength of the interaction between the slabs 
and the departure of the obtained frequencies from the results obtained with a single loop 
model. 



3.2. Resonant Damping of Oscillations 



Observations of transversa l coronal loop oscillations show that motions a re rapidly 
damped in only a few periods (lAschwanden et al.lll999l : iNakariakoy et al.lll999l). Corona l 
loops forming an arcade also show this property, such as reported by lVerwichte et al.l (120041 ). 
Resonant absorption o f wave energy ha s been proposed as a mechanism for the damping 
of loop oscillations by iHoUweg fc Yang (119881 ) . prior to observatio nal evidence, provided 
by TRACE observations, of damped transversal loop oscillations. iGoossens et al.l (120021 ) 
provided an interpretation to the damping observed in TRACE loop oscillations based on 
this mechanism. Although there is still no consensus on what causes the rapid damping 
of loop oscillations, several authors have presented comparisons between th eory and obser- 



vations that point out that resonant damping is a viable mechanism (see IGoossens et al. 
2OO2I : lAschwanden et al.l l2003l : lArregui et al.l l2007al . for example) . This damping mecha- 
nism relies on the coupling of global fast motions and localized Alfvenic motions due to the 
non-uniformity of the medium in the transversal direction. 



To study the damping of the normal modes described in § 13. 1| we now include non- 
uniform transitional layers that connect the internal and external densities at the edges of 
each slab {I 7^ 0). For the normal mode analysis of Equations ©-([Tj), a small but finite 
value of the resistivity has to be provided when computing the damping of oscillations. The 
magnetic Reynolds number in the solar corona is believed to be of the order of 10^^, hence 
the resistive diffusion time is very large. Such large values of Rm. are difficul t to h andle 



in any numerical scheme. On the other hand, as shown by iPoedts fc Kernerl (Il99ll ). the 
damping by resonant absorption becomes independent of resistivity, beyond a given value of 
Rm- In our numerical computations, resistivity has been chosen to be small enough for the 
imaginary part of the frequency to be independent of resistivity, so that we can be confident 
that solutions that correspond to the solar coronal large Reynolds numbers are obtained. 
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This condition has been checked to a high accuracy and for the solutions presented in this 
paper a value for the magnetic Reynolds number Rm = VAia/rj (with v^i the internal Alfven 
speed) in between 10^ and 10^ has been sufficient. 

We now consider the six different solutions presented in § 13.11 and compute their damp- 
ing due to the resonant coupling to Alfvenic motions. The same parameter values for the 
equilibrium configurations are considered, but now non-uniform transitional layers of thick- 
ness //a = 0.3 have been included. The results of our computations are displayed in Figured 
The left panels of this figure show the imaginary part of the frequency, for the six collective 
solutions under consideration, as a function of the position of the centers of the slabs. As 
can be seen, the splitting of solutions into symmetric and antisymmetric modes produces 
these modes to have a damping different from that of the corresponding single slab solutions. 
These differences go to zero in the limit of large separation between the slabs. As for the 
kink solutions, the antisymmetric mode has a stronger damping than its symmetric coun- 
terpart. This is also true for the sausage-like solutions. The most significant differences on 
the damping of the kink symmetric solution are produced for small distances between the 
slabs. However, the difference of the symmetric mode damping with respect to the kink mode 
damping of a single slab is below 10%. The kink antisymmetric solution damping, on the 
other hand, is more than double that of the kink mode of a single slab. Also significant are 
the differences between the sausage surface symmetric and antisymmetric damping and the 
corresponding sausage surface mode damping of a single slab. The imaginary parts for these 
two modes are 70% smaller and larger, respectively, than the value for the single slab mode. 
The right panels of Figure E] display the damping per period of the six solutions as a function 
of the perpendicular wavenumber, for a fixed value of the distance between the slabs. As 
we increase the value of perpendicular wavenumber the damping for both symmetric and 
antisymmetric solutions approaches the damping of the kink, sausage surface, and sausage 
body solutions of a single slab, in the limit of quasi-perpendicular propagation. Apart from 
the difference in the damping of these solutions with respect to those of the single slab and 
the dependence of the damping with perpendicular wavenumber, the most important result 
is that surface-like solutions, kink and sausage alike, are clearly strongly damped, within a 
few periods, for large values of ky, while sausage body solutions are unaffected by resonant 
couplings (see Figure O bottom row). The reason for the absence of resonant damping in 
the case of the two sausage body solutions is that their frequencies lie outside the Alfven 
continua produced by the presence of the four non-uniform layers. This leads to damping 
times, for these sausage solutions, that are of the order of the used Reynolds number and 
several orders of magnitude larger that the periods of oscillation. 

In order to better understand the results displayed in Figure [6] and to get physical 
insight into the reason(s) for the different damping rates of the different solutions, we have 
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analyzed the spatial distribution of the eigenfunctions for the six kind of solutions. Figures [71 
m and M display example eigenfunctions for the six types of solutions, for a particular set of 
parameters. The spatial distribution of eigenfunctions is rather similar to the ones shown in 
Figures [3|, IH and \5\ for undamped solutions, but now peaks at the non-uniform transitional 
layers are clearly visible in the perturbed transverse velocity of the kink and sausage surface 
modes, while they are absent in the case of the two sausage body modes. These peaks are 
an indication of resonant coupling to Alfven waves at the non-uniform layers. Aside from 
these peaks, the spatial distribution of the magnetic pressure perturbation has extrema at 
the four resonant layers in the case of the kink and the sausage surface modes, while it is 
maximum at the centers of the slabs in the case of the sausage body modes. The total 
pressure perturbation has no extrema at the resonant layers in the case of the sausage body 
solutions. The behavior of the eigenfunctions at the resonant layers and the funda r aenta l 
conservation laws that govern the resonant c ouplings were studied by ISakurai et al.l (Il99ll ) 
in ideal MHD and by iGoossens et al.l (119951 ) in dissipative MHD. These authors find that 
the derivative of the total pressure perturbation is zero at the resonant positions. This 
behavior is retrieved in the pressure perturbation for our damped solutions, but not in the 
two solutions for which resonant couplings are absent. 



4. SUMMARY AND CONCLUSIONS 



We have presented numerical results from a normal mode analysis for resonantly damped 
fast MHD oscillations in a system of two coronal loops. For simplicity, the adopted equi- 
librium model consists of two identical density enhancements in Cartesian geometry. The 
linearized resistive MHD wave equations are considered and solutions for obliquely propa- 
gating waves are obtained. The inclusion of oblique propagation together with non-uniform 
transitional layers that connect the internal and external densities produce the resonant 
damping of oscillations. 

The considered equilibrium supports a variety of different wave solutions. The normal 
modes of oscillation in a two-slab system are, in general, different to those of a single and 
isolated slab, due to the interaction between the slabs. The well-known fast kink and sausage 
normal modes, present in a single slab configuration, split into symmetric and antisymmetric 
solutions. We have concentrated our analysis on the six solutions coming from the splitting 
of the kink ar id sausage (body and surface) modes in a single slab with oblique propagation, 
described by lArregui et al.l (l2007d ). The frequencies of these six symmetric and antisym- 
metric solutions are in general different from the single slab s olutions. They depend upon 



the distance between the slabs (as shown by iLuna et al.l |2006| . for the kink symmetric and 
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antisymmetric solutions), but also on the magnitude of the perpendicular wavenumber. The 
inclusion of oblique propagation of perturbations has important effects on the dynamics of 
the system. First, it introduces two new surface-like solutions, with sausage parity, not 
present in the case of parallel propagation. Second, it decreases the effect of the interaction 
between the slabs, since solutions are more confined to the neighborhood of the density en- 
hancements. Both the distance between slabs and perpendicular propagation determine the 
difference in frequency between the single slab model solutions and the modes described in 
this paper. These differences are of the order of 15% at most and should, therefore, be taken 
into account in the modeling of these events. They are, however, well below the observational 
uncertainties for t he measured periods in transversal coronal loop oscillations, which are of 
the order of 40% (lAschwanden et al.ll2002l ). 

Next, we have studied the damping of the different solutions due to resonant coupling 
to Alfvenic motions. The period and damping rate of the six solutions have been computed 
and the spatial structure of eigenfunctions analyzed in order to find some physical insight 
into the obtained results. We find that surface-like solutions, both kink and sausage, are 
strongly damped. On the other hand, the two sausage body solutions are unaffected by 
resonant absorption. This is due to the absence of resonance due to the fact that the 
frequencies of these solutions are outside the corresponding Alfven continua. A comparison 
with the single slab model results reveals that differences in damping are important for 
small values of the distance between the slabs. The most significant differences arise for the 
kink antisymmetric solution and the two sausage surface solutions, with differences in the 
imaginary part of the frequency of the order of 125% and 70% with respect to the damping 
of the corresponding single slab modes. These differences are of the order or larger than 
the observational uncertainties for t he measured damping ti mes in transversal coronal loop 
oscillations, which are around 60% (lAschwanden et al.ll2002l ). 

In this work a magnetic two-slab system has been used to model a two-loop structure. 
Some of the results found in slab geometry cannot be translated to cylindrical coronal loop 
systems, while some others give a good approximation to the phase speeds in cylinders. For 
instance, the sausage surface solutions described in this paper are not present in a system 
of two cylinders. Also, the sausage body solutions would correspond to perturbations in 
each cylinder with m = 0, thus they cannot be resonantly damped solutions. On the other 
hand, slabs are known to be poor wave-guides in comparison with cylinders. However, our 
results indicate that the inclusion of oblique propagation produces a sharper drop-off rate of 
the eigenfunctions in the external medium, resulting in a confinement similar to the case of 
cylinders and a good approximation of the phase-speed of the kink mode in a slender flux 
tube. The normal mo des of oscillation i n a system of two homogeneous cylinders have been 
computed recently by iLuna et al.l (120071 ). These authors find four solutions that correspond 
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to transversal oscillations, symmetric and antisymmetric, in the plane containing the two 
cylinders. The frequency of these solutions is seen also to differ from the frequency of the 
kink mode for a single cylinder for small distances between the cylin ders. A natural ext ension 



to the study presented in this paper and to the results described by iLuna et al.l (120071 ) is the 
analysis of the damping in a similar two-cylinder configuration. 
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Fig. 1. — Schematic representation of the two density enhancements (hght-shaded regions) 
of half-width a, centered at ±xo representing a system of two coronal slabs in the direction 
transverse to the equilibrium magnetic field. These enhancements with internal density pi 
connect to the external medium, with density p^., by transitional non-uniform layers (shaded 
regions) of thickness /. 
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Fig. 2. — Left, Frequency as a function of the position of the slabs centers for different values 
of the perpendicular wavenumber for the kink solutions (top), the sausage surface solutions 
(middle), and the sausage body solutions {bottom). Right, Frequency as a function of the 
perpendicular wavenumber for the same six types of solution, for a fixed distance between 
the slabs, Xq = 5a. In all the figures the solid lines correspond to the symmetric solution 
and the dashed lines to the antisymmetric solution. 
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Fig. 3. — Real part of the transversal velocity component, Vx, and imaginary part of the 
perturbed total pressure, Pt, for the fundamental symmetric (left) and antisymmetric (right) 
kink solutions, for several values of the perpendicular wave number: kyU — (dashed lines), 
kytt = 0.5 (dotted lines), and kya = 1 (solid lines). Note that the character of the eigen- 
functions for the symmetric mode changes from body-like to surface-like solution as ky is 
increased. In all figures the light-shaded regions represent the density enhancements. These 
solutions have been obtained with a uniform computational grid with = 10 000 points 
in the range —80 < x/a < 80 (except for the kink solutions with ky = for which a range 
-200 <x/a< 200 has been used). 
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Fig. 4. — Same as Figure [3] for the sausage surface, symmetric [left) and antisymmetric 
[right) solutions and two values of the perpendicular wavenumber: kya = 0.3 (dotted lines) 
and kya = 0.5 (solid lines). These solutions are surface-like for all values of the perpendicular 
wavenumber. 
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Fig. 5. — Same as Figures [3] and H] for the sausage body, symmetric (left) and antisymmetric 
(right) solutions and two values of the perpendicular wavenumber: kya = 0.5 (dotted lines) 
and kytt = 1 (solid lines). Note that for these solutions Vx is very similar to the ones depicted 
in Figure HI while the pressure perturbation in and around the slabs differ substantially. 
These solutions are body-like for all values of the perpendicular wavenumber. 
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Fig. 6. — Ze/t, Imaginary part of the frequency as a function of the position of the slabs 
centers for the kink solutions {top)^ the sausage surface solutions {middle)^ and the sausage 
body solutions {bottom) for kya = 0.3 and l/a = 0.3. Right, Damping time over period as a 
function of the perpendicular wavenumber for the same six types of solution, for the same 
parameter values and a fixed distance between the slabs, Xq — 5a. In all the figures the 
solid lines correspond to the symmetric solution and the dashed lines to the antisymmetric 
solution. A value for the magnetic Reynolds number Rm — VAia/f] — 10^ has been used in a 
non-uniform grid with A^^, = 14 000 points (2 000 in each of the four non-uniform transitional 
layers) in the range —80 < x/a < 80. 
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Fig. 7. — Real part of the transversal velocity component, f^,, and imaginary part of the 
perturbed total pressure, Pt, for the fundamental symmetric (left) and antisymmetric (right) 
kink solutions, for kytt = 0.3 and non-uniform transitional layers with thickness l/a = 0.3. A 
value for the magnetic Reynolds number Rm = VAia/f] — 10^ has been used in a non-uniform 
grid with — 14 000 points (2 000 on each of the four non-uniform transitional layers) 
in the range —80 < x/a < SO. In all figures the hght-shaded regions represent the density 
enhancements and the shaded regions the non-uniform layers. 




Fig. 8. — Same as Figure [7] for the fundamental symmetric {left) and antisymmetric {right) 
sausage surface solutions. 
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Fig. 9. — Same as Figures [7] and [8] for the fundamental symmetric (left) and antisymmetric 
(right) ge body solutions. 



